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In the context of the imaginary-time formalism for a scalar thermal field theory, it is shown 
that the result of performing the sums over Matsubara frequencies associated with loop Feynman 
diagrams can be written, for some classes of diagrams, in terms of the action of a simple linear 
operator on the corresponding energy integrals of the Euclidean theory at T — 0. In its simplest 
form the referred operator depends only on the number of internal propagators of the graph. 

More precisely, it is shown explicitly that this thermal operator representation holds for two 
generic classes of diagrams, namely, the two- vertex diagram with an arbitrary number of internal 
propagators, and the one-loop diagram with an arbitrary number of vertices. 

The validity of the thermal operator representation for diagrams of more complicated topologies 
remains an open problem. Its correctness is shown to be equivalent to the correctness of some 
' diagrammatic rules proposed a few years ago. 

(N ' 
O 

Q 

I. INTRODUCTION 

In the imaginary-time formalism, the calculation of a loop diagram in quantum field theory at finite temperature 
necessarily involves sums over Matsubara frequencies 0, an operation that we shall generically call the Matsubara sum 
associated with the graph. Although this sum can be computed in a number of ways, usually in a systematic fashion, 
such computations can become quite tedious for higher loop diagrams |2t |3j. 

In reference a set of simple diagrammatic rules were postulated to write down an explicit expression for the 
result of performing the Matsubara sum associated with any finite-temperature Euclidean Feynman graph (in a scalar 
C*~) 1 theory). Because of the similitude of the diagrammatic expansion with the one associated with the non-covariant 
old-fashioned perturbation theory formalism (at zero temperature), these diagrammatic rules will be referred to as 
the OFPT-rules. 

i' Although in reference Q the OFPT-rules were explicitly verified to hold for a few nontrivial diagrams, they were 
' , presented as a sort of empirical discovery, with no rigorous proof given. 

In this paper we restate the diagrammatic analysis of reference [4J in an algebraic rather than diagrammatic fashion, 
and extend its validity to two particular classes of diagrams, to be described below. For these diagrams we establish 
that the full result of performing the Matsubara sum associated with a given Feynman graph can be completely 
determined from its zero-temperature counterpart, by means of a simple linear operator, as shown in below. We 
have termed this result the thermal operator representation (TOR) of the Matsubara sum. 

The two classes of diagrams for which we have been able to prove the correctness of the TOR are: (a) diagrams 
with two vertices and an arbitrary number of scalar internal propagators; and (b) one-loop diagrams with / vertices 
and / scalar internal propagators, with I > 1. In what follows, whenever we refer to a Feynman diagram we implicitly 
assume that the diagram actually belongs to one of the classes just described, except when specifically qualified 
otherwise. 

The precise mathematical formulation of the thermal operator representation is presented in the next section. 
Leaving out many of the technicalities, its contents is as follows. Consider the Matsubara sum of a (amputated) 
scalar loop Feynman graph with / internal lines and external Euclidean 4- momenta P a = (p Q ,p a ). [A word about 
the notation: in order to avoid clutter, we will omit the customary superscript on Euclidean energy variables. Since 
we shall not denote in this paper the modulus of a 3-momentum vector p with the corresponding italic symbol p, 
there should be no danger of confusion]. Instead of following the usual practice of parameterizing all internal line 
4-momenta in terms of a few independent loop 4-momenta, by explicitly requiring 4-momentum conservation at each 
vertex, we choose to assign each internal line independent 3-momentum and Matsubara frequency fcj and impose 
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4-momentum conservation by means of an appropriate number of delta functions. In this form, the Matsubara sum 
will depend only on the external Euclidean energies p a (which enter through Kronecker delta functions enforcing 
energy conservation at each vertex), on the kinematic energies of the internal lines, := (k| + mf) 1 ^ 2 appearing in 
the propagators and, of course, the temperature T. Since there is no explicit dependence of the Matsubara sum on 
spatial 3-momenta, external or internal, we shall suppress all reference to these in this paper, whenever possible. 

Let the unsubscripted symbols p and E denote, respectively, the full set of Euclidean external and kinematic 
internal energies, p := {p\,p2, ■ ■ ■ ,p n } and E := {E\, E%, . . . . E{\. Now, if we introduce the Matsubara D-function of 
the graph, D(p, E, T), essentially as the Matsubara sum multiplied by the product of all internal kinematic energies, 
then we claim that 

D(p,E,T) = 6(E,T)D (p,E), (1) 

where 0(E,T) (whose explicit form we give in the next section) is a linear operator that depends on the topology 
of the diagram but is independent of the external Euclidean energies p. The object acted upon by this operator, 
Dq(p,E), is simply the corresponding D-function for the Euclidean zero-temperature graph, Dq(lv,E), defined for 
real and continuous external energies uj :— {u>\, uj 2 , . . . , w n }, evaluated at uj = p: 

D Q (p,E) = D {uj,E) . (2) 



We shall call 0(E,T) the (Euclidean) thermal operator. 

As we shall see in the next section, the thermal operator has a form that can be readily and naturally extrapolated 
to diagrams of arbitrary topologies. Although this makes it tempting to conjecture that the thermal operator repre- 
sentation holds for completely arbitrary diagrams, this remains an open problem and more work is needed to settle 
the issue. 

However, if true in general, the representation Q would have several immediate important consequences: (a) it 
would show that the full finite temperature result is encoded in the zero-temperature function D (p, E), rendering the 
actual computations of the Matsubara sums completely unnecessary; (b) since all dependence on external energies p 
is contained into the zero-temperature function Dq, any analytic continuation of D{p 1 E, T) to complex values of the 
external energies, physically meaningful or not, would need only be carried out on Dq. By the same token, the study 
of imaginary parts of analytically continued Euclidean Green functions, i.e., the subject of cutting rules, would need 
only be done at the level of the zero-temperature function Do, since the thermal operator is real (we give an example 
of this in the last section of this paper), (c) since the thermal operator is bounded as the internal energies E\ tend 
to infinity, it would be enough to renormalize Dq in order to renormalize the full finite temperature result. This is 
consistent with a well-known result in renormalization of thermal field theories. 

Although there have appeared in the literature several works that touch upon the relationship between the full 
calculation of finite-temperature Feynman graphs and their zero-temperature counterparts (usually interpreted in 
terms of forward scattering amplitudes in vacuum), both in the Euclidean imaginary-time [5( and in the real-time 
formalisms@, we are unaware of any discussion of a representation of the simple form as given here. 

We emphasize that all the results presented in this paper are formulated in the context of the Euclidean imaginary- 
time formalism, and we will have nothing to say here about their relationship to or consequences for the real-time 
formalism, except for the remark made above about the possible analytic continuations of the Euclidean Green 
functions to complex values of the Euclidean external energies. The latter subject has been studied at great length in 
the literature 0, along with the connection between different analytically continued Euclidean Green functions and 
the retarded, advanced or time-ordered Green functions of the real-time formalism, and the subject of cutting rules 
in the real-time formalism^. 

The structure of the paper is as follows: In section [D] we shall present the general form of the thermal operator 
representation (TOR) for the Matsubara sum of a general scalar graph, in two alternative forms. In sections lTTTl HVI and 
[V]we prove that the TOR holds, respectively, for a one-loop single-propagator tadpole-like graph, for a generic graph 
with two vertices, and for a generic one-loop graph; the number of internal propagators is allowed to be arbitrary (but 
at least equal to two) in the last two cases. Additional supporting evidence for the validity of the TOR for graphs of 
arbitrary topologies and our conclusions are given in section IVll The reformulation of the old-fashioned perturbation 
theory rules of reference jj] in the form of the present representation has been relegated to an appendix. 
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II. A REPRESENTATION FOR THE MATSUBARA SUM 



In a scalar field theory, the mathematical expression corresponding to an amputated graph with n + 1 vertices 
(n > 1), I internal lines, and external 4-momenta P a — {jp a , p a ) has the form 



i: " Al +1 / [ n?^S^ n^w^«,n (3) 

i—i \ ) 1 y— i 
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where A represents the coupling constant and S is the symmetry factor of the graph; is the spatial 3-momentum 
of the i-th internal line and E{ := (kf + mf) 1 ^ 2 is its associated kinematic energy; kv denotes the total 3-momentum 
entering vertex V; the unsubscripted symbols p and E denote, respectively, the full set of Euclidean external and 
kinematic internal energies, p :— {pi,P2, ■ ■ ■ ,Pn} and E := {E\, E<i, . . . , E{\; and T is the temperature. The delta 
functions ensure conservation of spatial 3-momentum at each vertex, so that the integration measure reduces essentially 
to an integration over the 3-momenta of the L = I—n independent loops. In the finite temperature Euclidean formalism 
all lines, external and internal, carry discrete Euclidean energies which are integer multiples of 2itT. Each internal 
line has an associated Matsubara frequency, denoted by kj = 2irTnj. The D-function is given by the normalized 
Matsubara sum 

i 

D{p,E,T)= lE T L J2 H A ( k 3> E MP,k), (4) 

n 1 ,n 2 ,...,n I j=l 

where 

lE ~\[2E l . i (5) 

i=l 

L is the number of independent loops in the graph, and A(kj,Ej) is the scalar propagator associated with the j-th 
internal line, with 

A(*,£):=^. (6) 

The sums over each rij run from — oo to +oo. The (^-function, with k = {k\, . . . , kj}, is a generalized Kronecker 
delta which ensures conservation of energy at each vertex. The topology of the diagram is totally contained in this 
generalized delta. 

The OFPT-rules given in which are reproduced in appendix IVIII Al were conjectured to allow us to write down 
the complete result for Q by a simple diagrammatic analysis. But as shown in appcndix lVIII til there exists a simple 
algebraic representation for the diagrammatic OFPT-rules, so that the conjecture of reference Q can be recast in the 
following terms: 

Statement 1. [Thermal Operator Representation] The D-function defined in (@J for an amputated Feynman 
graph can be expressed in the form 

D(p,E,T) = 6(E,T)D (lo,E) , (7) 

LJ—P 

where Dq(u),E) is the D-function of the Euclidean zero-temperature graph and 0(E,T), the thermal operator, is the 
following linear operator: 

I 

6(E,T) := l+^2m(l+Si) + ^ n n n l2 (l +5 n )(l+5 l2 ) 

l in ' l2) (8) 

+•••+ e' n^( i+ ^)- 

1=1 

Here ni = n(Ei), where n(E) — (e^ E — l) is the Bose-Einstein thermal occupation factor: Si :— S Ei is a reflection 
operator, S x f{x) := /(— x); the indices i\, i2, ■ ■ ■ run from 1 to I (the number of internal propagators) and the symbol 
. . . ,ik) stands for an unordered k-tuple with no repeated indices, representing a particular set of internal lines. 
The primes on the summation symbols imply that certain tuples (ii, . . . ,ik) are to be excluded from the sums: those 
such that if we snip all the lines i\, . . . , iu then the graph becomes disconnected. 



FIG. 1: A one-loop single-propagator diagram. 



Note that the operator 0(E, T) contains products of at most L thermal occupation factors n(Ei), since for a L-loop 
graph the maximum number of lines that can be snipped without disconnecting the graph is precisely L. This generic 
feature of the thermal graph in the imaginary-time formalism is of course well known. However, as discussed in 
sections llVI and Ivl there exists a simpler algebraic form for the thermal operator: 

Statement 2. [Simpler form of the Thermal Operator] When acting on the zero-temperature D-function, 
D$(p,E), the thermal operator 0(E,T) can be replaced by the the simpler 

i 

6+{E,T)=Y[[l + n l (l + S l )}. (9) 

i=l 

Note that the operator O* (E, T) in @ can be expanded as in iJSJ , except that the summation symbols carry no 
primes, that is, all tuples (ii, . . . , ik) (1 < k < I) are allowed in the sum. Clearly, the form © will follow from ||SJ| if we 
can somehow show that tuples associated with disconnected graphs (the ones excluded from the summations in (JSJ) 
give rise to operators that produce a vanishing contribution to the D-function in Q. So, the simpler representation 
will follow from JSJ if the following statement is true: 

Statement 3. [Cut sets do not contribute] The zero-temperature D-function, Dq(lo,E), is annihilated by the 
operators 

A{C):= 11(1 (10) 

hec 

where C stands for a cut set of the graph, that is, any set of indices i\, . . . , ik such that the graph becomes disconnected 
if the corresponding lines are snipped. 

We make clear at this point that, although we make reference to cut sets, we imply no connection to the concepts of 
cuts and cut diagrams as they are usually understood in diagrammatic quantum held theory. Cut sets are determined 
solely by the topology of the diagram, and have no further mathematical or physical meaning. 

The goal of the next three sections is to prove that these statements are indeed true for the two generic types of 
graphs described in the introduction. The strategy of the proof will be to evaluate the Matsubara sums contained in 
D(p, E, T) by conventional means, namely the contour integration method or the Saclay method, and then show that 
the result can be written as in the right-hand side of Q. 



III. THE SIMPLEST LOOP DIAGRAM 

We begin by considering a one-loop graph with only one internal propagator, as the one shown in hgurc ^ This 
particular graph contributes at first order to the self-energy in the A</> 4 theory. The actual number of external legs of 
the graph is unimportant, since we are only interested in the Matsubara sum associated with the loop. Although we 
could have considered this graph as the simplest case of the generic one- loop graph considered in section [V] we prefer 
to analyze it separately, since the proof given in section applies more naturally to the case of two or more internal 
propagators. 
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According to the ^-function for the graph of figure is simply given by 



D( P ,E,T) = (2E)T f {27rTn ] 2 + E2 - (ID 



n— — oo 



The sum above can be computed in a variety of ways and the result is well known lj . One obtains 

D(p,E,T) = l + 2n(E), (12) 

where n(E) — (e^ E — l) 1 is the Bose-Einstein thermal occupation factor. The zero-temperature D-function, 
Do(uj,E) can be computed directly from its definition, 

oo 

Do( ., E ) = <2 E >/§£^, (13) 

— OO 

or simply by taking the limit T — > of D(p, E, T) in l(T2|). The result is 

D (u,E) = l. (14) 

Since a constant function is unchanged by the reflection operator S E defined by 

S E f(E) := f(-E), (15) 

where / is any regular function in the variable E, we certainly have the identity 

D( P ,E,T)= [l + n(E)(l+S E )]D (p,E), (16) 

which proves that the thermal operator representation given by (0) and (jSJ does hold for the simple graph we are 
considering. 

IV. THE TWO- VERTEX DIAGRAM 
A. Calculation 

The Matsubara sum for the two- vertex diagram with / internal propagators shown in figure El is most conveniently 
calculated using the Saclay method|2|, which we now briefly review. 

Let K :— (fc,k) be the Euclidean 4-momentum vector associated with a given internal line; k is a Matsubara 
frequency to be summed over. 

Then each scalar propagator, 

A <*> := k*t^ = wh? : = A(fc <^ (17) 



where E :— E^ = Vk 2 + m 2 , is represented as 



A(k,E)= [ dTe lkT A(r,E), (18) 
Jo 

where (3 = 1/T, as usual. The mixed propagator A.(t,E), < r < (3, is given by 

A(r, E) = -L [ (1 + n{E)) e ~ ET + < E > ET ] > ( 19 ) 

where n(E) = {eP E - l) For our purposes, it will be convenient to use the following representations for the mixed 
propagator ljT51) : 

A(r,£) = ±[l + n(E)(l+S E )]e- E \ (20) 

= ^^[l + e-^ E S E ]e-^, (21) 



FIG. 2: A two-vertex diagram with / internal lines. 



where Se is the reflection operator defined in l|15|l . Substituting the representation (|2U|) back into (|18|l and using the 
fact that the operator (1/2E)[1 + n(E)(l +Se)] is linear, we obtain the following equivalent Saclay representation for 
the scalar propagator: 

A(fc, E) = ± [1 + n{E )(l + S E )} J P dr e^~ E ^ . (22) 



Consider now the two-vertex diagram with I internal lines of figure Let P = (p, p) its external (incoming) 
4- momentum (note that here p stands for a single Euclidean energy variable) and let Kj = (kj = 27rTn :) , kj), j — 
1, . . . ,1, be the 4-momenta of the internal lines, flowing from the left to the right vertex. The Matsubara D-function 
corresponding to this graph is given by 



D^E^^eT 1 - 1 J2 \{^{k 1 ,E 3 )5 kl+ ... +kl , p , 

ni,Ti2,...,ni j = l 



(23) 



where the delta function is a Kronecker delta enforcing conservation of energy at both vertices, Sj=i kj — P, and je 
was defined in (JSJ). 

Now, because the variables p and kj are quantized in units of 2ttT, the Kronecker delta in 12.'il) can be represented 

as 



6 kl+ ... +kl , p = T / dre -i-(fci+-+fex-P) ) 



(24) 



so that the sums over the integers rij decouple: 



D(p,E,T)= lE T I f 'dre^TT 



—irkj 



(25) 



Using now the Saclay representation (122[) for each propagator A(kj,Ej) (with integration variable tj) we find 



D(p, E, T) = T 1 [] [1 + n(Ej)(l + S E] ] 
r0 



x 



/ dre ipT TT / dT 3 e~ ET i Ve'' T '- r »*' 
J ° ,=i J o 



(26) 



But 



rij n 

= 5(rj - r), for < Tj.r < f3. 



(27) 
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so that the final result for the Matsubara D- function for the graph of figure [5] is: 

' f/3 



D(p, E,T) = TT [1 + n(S,-)(l + 5/)] / dre^-***) 
„■_-, Jo 



3=1 
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= J] [1 + n(^-)(l + 5,)] . _ p , (28) 

j—l °P -^tot 

where <Sj := Se,, E to t ■— Y2j=i Ej and we have used the fact that exp(i/3p) = 1. 

B. Proof of the Thermal Operator Representation 

We will now show that the result l|28fl can be put into the form Q , where the zero-temperature D-function for the 
graph of figure [21 is given by 



D (p,E) = - 



1 



1 



ip - E tot ip + E t , 



(29) 



as can be easily be obtained from a calculation in the old-fashioned perturbation theory formalism. First, we observe 
that this function satisfies statement 3. In fact, since the only cut set of the two- vertex diagram is the set of all lines, 
we only need to show that the function is annihilated by the operator 



A:=H(1+S 3 ). 

3=1 



(30) 



But since S x is a reflection operator (S% = 1) we have 



3=1 



1 



ip - E t , 



3=1 
I 



3 = 1 



ip + Etot ' 



(31) 



so that indeed ADoip, E) = 0. Therefore it is enough to show that Q holds with the thermal operator in the form 
@. But this follows immediately from the following identity: 



3 = 1 



-0Eu 



ip - E t , 



-0E t , 



ip - E u 



n^[i+ g -^] 

3=1 
3=1 

J 

IJ[l+»(Si)(l + 5 J )] 



1 



3 = 1 



1 



where we have used the property (e ^ El Sj) = 1- 



V. THE ONE-LOOP DIAGRAM 



A. Calculation 



The calculation of the Matsubara sum for the one-loop diagram with / vertices and / internal propagators shown 
in figure |2| is most conveniently done using the standard contour integration method^. If a meromorphic function / 
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has no singularities along the imaginary axis and k stands for the Matsubara frequency k = 2nTn, then 

T fW =2-9 tt z )< z ) dz ' ( 32 ) 

n=-oo •' C 

where n(z) = (e^ z — l) and C is the positive contour that runs vertically upwards the complex z-plane, infinitesi- 
mally to the right of the imaginary axis, from e — ioo to e + ioo, and then comes back vertically and infinitesimally 
to the left of the imaginary axis, from — e + ioo to — e — ioo, with s = + . 

If \n(z)f(z)\ goes fast enough to zero when \z\ goes to infinity we can change the contour of integration to two 
negatively oriented semicircumferences, one on each side of the imaginary axis, with radii tending to infinity. Thus, 
by Cauchy's integral theorem, 



t J2 /(*"*) = -E Res *=*» [/(*m*)]. 



(33) 



where zi are the poles of the function f(z). 

Consider now the one loop graph of figure |3 Let Pi — (pi,Pi) be the external incoming momenta at each vertex. 
Letting k = 2imT be the Matsubara frequency of line 1, the Matsubara D-function in this case can be reduced to the 
form 



D(p, E, T) = 1E TJ2 A(fc, E 1 )A(k + pa, E 2 ) ■ ■ ■ A(k + J^Ph^i), 

n j=2 



(34) 



where the energies Ei are defined as before. Introducing a new set of variables Uj := X)/=i Pi ~ Pi (j = 2, . . . , I) and 
letting u\ := 0, we can write as 



D{p, E, T) = 1E TJ2 II A ( fc + u i' E i)- 

n j = l 



Next using the identity 



we get 



k 2 + E 2 2E \ik + E ik-E 



2E ^ ik-o-E' 



(35) 



(36) 



D(p,E,T)=(-iyTj2\l[ 

n I j=l 



(37) 
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where now a :— {a±, 02, ■ ■ ■ , 07}. If the function between brackets in l|37|) is called f(ik), then we see that the poles 
of f{z) are located at z\ = — iu\ + &iEi, so that the application of ij-i-il) gives us 



D(p,E,T) = (-1)™ J2J2 n ^ E ^E 



a 1 = 1 



- ui) + aiEi - UjEj 



(38) 



In order to express the result for the /^-function in terms of Bose-Einstein factors of positive argument only, we will 
perform the summation over ai explicitly. Introducing the notation <jj :— (a±, . . . , <tj_i, ■ ■ ■ , <Ji) and using the 
identity n{—E{) = — (1 + n(E t )) we find 



d(p^t)=(-d^j:e n? 



+n[E l ) 



i(uj - ui) -Ei- <7 3 Ej 



n 



n 



^ i(uj - ui) -Ei- ajEj i( Uj - ui) +Ei- ajEj 



(39) 



In terms of the auxiliary function 



di(p,E) =J2U 

=nz 



i(uj - ui) - Ei - UjEj 



a _2 

— ui) — Ei — UjEj 



(40) 



and the reflection operator, Si :— defined in l|15l) we have 



D(p,E,T) = [di(p,E) +n(Ei)(l + S E ,)d,(p, E)] 



(41) 



1=1 



B. Proof of the Thermal Operator Representation 

We shall prove now that the D-function ljl"T)l for the one-loop graph of figure can be written in the form Q, as 



D(p,E,T) 



l + ^n(^)(l+5 J 

3=1 



D (p,E), 



(42) 



with 



D (p,E) = (-l) I + 1 J2di(p,E). 



(43) 



1=1 



Since the graph of figure gets disconnected if two or more lines are snipped, the thermal operator has terms no 
higher than linear in the Bose-Einstein factors n(E). But equation l|42l) will reduce to eqn. H41|l if the operator (1+iSj) 
annihilates the auxiliary function di(p, E) when j ^ I. This is indeed the case: from equation (|40(l we see that when 
3+1 



d«(p,£0=E 



1± TT \ "* a J_ 

i(uj - ui) - Ei - ajEj if .-^ i(uk ~ ui) ~ E - a k E k '' 



nz- 



(44) 



i(uj - ui) - Ei + ajEj Af i(uk - uA — Ei — a k E k 

' J J k^=l] (Jfc v ' 



Sjdiip, E), 
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which means that 

(l + S j )d l (p,E) = 0, i£j?l. (45) 

Statement 1 is then valid for the one- loop graph of figure El 

Furthermore, statement 3 is also true for this graph. In fact, any cut set will at least contain two lines, say lines i 
and j. But then 

(l + S i )(l+S j )d l {p,E)=0, (46) 

since, for any given /, either i or j will be different from I (since i ^ j), leading to a vanishing contribution because 
of gjjjl. 



VI. FURTHER EVIDENCE AND CONCLUSIONS 



One piece of evidence in favor of the general validity of the representation {JJ is provided by a comparison with a well 
known result of thermal field theory, first formulated by Weldon[lfT|. concerning the interpretation of the imaginary 
part of the retarded self-energy Tlji in terms of the direct and inverse decay rates of a particle propagating in the 
thermal medium. A well-known result of quantum statistical mechanics |lj is that the full retarded self-energy Hr 
can be obtained from the Euclidean self-energy by analytic continuation as 



n fl (w, p) = -Ilf)(i(u + ie), p), 



(47) 



where u> stands for a real continuous variable. In the context of perturbative quantum field theory, the imaginary part 
of 11^ is given in the form of integrals over phase space of amplitudes squared, weighted by certain statistical factors 
that account for the possibility of particle absorption from the medium or particle emission into the medium^j- F° r 
example, for the one- loop 2-vertex diagram corresponding to figure 0] in the appendix that follows, the result for the 
imaginary part of the retarded self- energy is (we have set g = 1) 



ImIIfl(cj,p) = 



d 3 k 



1 



(2tt) 3 AE1E2 



(1 + m + n 2 )(S{uj — Ei — E 2 ) - 6(ui + E X + E 2 )) 
-(ni - n 2 )(5(Lu -E x + E 2 ) - 5{lo + E x - E 2 )) 



(48) 



where n, = n(-Ej). But from the general form ijJJ for a diagram in the Euclidean formalism, it is clear that the 
imaginary part of the analytically continued diagram is determined by the analytic continuation of its D-function. 
The general validity of our main representation in the form (0 would imply that the latter is in turn completely 
determined in terms of the analytic continuation of the zero-temperature D-function, Dq, since the thermal operator 
O is real and does not involve the external momenta. 

For the particular simple diagram we are considering, which is actually a special case of the general 2-vertex graph 
considered in section llVl the Thermal Operator Representation has been proven to hold. Hence 



Im D + ie) , E lt E 2 ,T) = 6 (E u E 2 , T) Im D (i (w + te),E 1 ,E 2 ). 



(49) 



The last imaginary part could in principle be obtained from the standard cutting rules that apply in zero-temperature 
field theory, without having to compute _D itself. In this case, however, we have the closed result (|29f) for _D , which 
allows us to compute directly 



ImA) (i(uj + ie),E 1 ,E 2 



Im 



1 



1 



uj + Ei + E 2 + ie uj — Ei — E 2 + ie 
= -7T [S(u + Ei+ E 2 ) - 5(lu -Ei- E 2 )] . 



Now in this case the thermal operator is given by 
6{Ei,E 2 ,T) = l + ni(l+5i) 

Since 



n 2 {\ + S 2 ) = 1 + m + n 2 + niSi + n 2 S 2 . 



(50) 



(51) 



mSi [S(u + Ei+ E 2 ) - S(ui -Ei- E 2 )] = m [5(u - E x + E 2 ) - 8{uo +E X - E 2 )} , 



11 



etc., we readily obtain 

OlmDo = TTUl+ni+n 2 )(6(u-E l -E 2 )-6(u + E 1 +E 2 )) 

-(»i - n 2 )(6(LU -Ei + E 2 ) - 5{uj + E x - E 2 ))] , (52) 

thereby reproducing l|48(l . with all the correct signs and thermal factors. 

In this paper we have restricted our attention to some simple diagrams in the finite-temperature imaginary-time 
formalism for a scalar relativistic field theory. We have shown that the full result of performing the Matsubara sum 
associated to any given Feynman graph can be obtained from its zero-temperature counterpart, by means of a simple 
linear operator. Given the general form JSJl of the thermal operator, which can be readily and naturally extrapolated 
to diagrams of arbitrary topologies, it is not at all implausible that the representation Q be actually valid in complete 
generality. This generalization remains an open problem however, and work in this direction is in progress. 

An analysis similar to the one presented here should apply in a theory containing fermions; the algebra will be 
slightly more complicated because of the spin structure. We have deferred this analysis, as well as the extension of 
our results to gauge theories, until we have been able to prove or disprove that the Thermal Operator Representation 
put forward in this paper does indeed hold for an arbitrary loop graph in a scalar field theory. 
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VIII. APPENDIX 
A. The OFPT-rules 

The rules originally put forward in reference £| to write down an explicit expression for the Matsubara D-function 
corresponding to the general scalar graph considered in section [n] are given by the following statements (it might be 
useful to refer to figure 0] at this point): 

a. For each external line, characterized by a real Euclidean 4- vector (pi,pi), define its energy as ipi. For each 
internal line define its energy as Ei = (k^ + mf) 1 / 2 , where is the 3-momentum carried by the line and m, is 
the mass of the propagating particle. 

b. Define a direction of time or energy flow (which we shall take conventionally from left to right) and consider all 
possible orderings of the vertices along this direction (see, e.g., Figs.0]a and^Jb. For a graph with n + 1 vertices 
there will be (n + 1)! such orderings). 

c. For each time-ordered graph generated in (b) consider, in addition to itself, all possible connected graphs that 
can be obtained by snipping any number of internal lines. Each line that is snipped becomes a pair of legs 
we shall call thermal legs. Attach a cross to their ends to distinguish them from the original external lines of 
the graph. Both legs of a given pair inherit the energy Ei of the internal line that originated them. However, 
one leg must be oriented as incoming with energy Ei and the other as outgoing with energy Ei. Both possible 
orientations have to be considered, each one generating a different diagram (see, e.g., Figs. 0]c and^Jd). 

d. For each graph in (c), define its total incoming energy, Ei nc , as the sum of all incoming external energies plus 
the energies of all incoming thermal legs that join the diagram before their outgoing partner (e.g. as in Figs. 
0]d and0]f). Thermal leg pairs that satisfy this property shall be referred to as external, and those that do not 
as internal (e.g. as in Figs. 0Jc and0]c). Then, associate to this graph an expression equal to the product of 
the following factors: 
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draw a full vertical division (a "cut" ) between each pair of consecutive time-ordered vertices (there are n 
such cuts in a graph with n + 1 vertices); for each cut include a factor 

(53) 



Ein.r E CA < 



where E cut is the total energy of the intermediate state associated with the cut, defined as the sum of the 
energies of all the lines that cross the cut in question (as in zero-temperature time-ordered perturbation 
theory), plus the energies of all internal thermal pairs whose originating internal line would have crossed 
the cut. 

2. include a thermal occupation factor = n(Ei) for each thermal pair (of energy Ei) in the diagram (if 
any). 

3. include an overall factor of (—1)", where n + 1 is the number of vertices. 

e. The integrand D(p,E,T) in i.e., the Matsubara D-function of the graph, is the sum of the expressions 
computed according to rule (d), over all the graphs in (c). 



B. An algebraic approach to the OFPT-rules 

Let us call D^(p, E,T) the expression for the D-function generated according to the OFPT-rules. A trivial check 
the OFPT-rules do satisfy is that they yield the known correct result in the limit T — > 0, keeping the external 
Euclidean energies p fixed. In fact, in the limit T — > all the thermal factors n(E) vanish, so that according to the 
rules Dr(p,E,0) is just given by all possible time ordered diagrams with no snipped lines, calculated according to 
rule (d) above. But this is precisely the result one would obtain calculating the T — Euclidean graph (with external 
momenta (pupi)) using old-fashioned perturbation theory |9j. We have therefore: 

D o (p,E) = D R (p,E,0), (54) 

where D$(p,E) is the D- function associated with the zero-temperature Euclidean Feynman graph. Hence, the rules 
hold at T = 0. 

At finite temperature, we get extra contributions according to rules (c) and (d) above. Now, instead of considering, 
as commanded by rule (c), all possible connected graphs that can be obtained by snipping any number of internal 
lines of a given "un-snipped" time-ordered graph, let us rather group the snipped diagrams according to which lines 
are snipped, regardless of the time-ordering. Take for instance all the diagrams which have only the z-th line snipped 
(i is fixed). A set of this type is conformed, for instance, by the diagrams (c) to (f) of Fig0J It follows directly 
from rule (d.l) that a diagram in which the snipped line forms an internal thermal leg pair (i.e., we have a "closed" 
snipping, as in Figs. 0]c and 0Je) has exactly the same mathematical weight as the zero-temperature "un-snipped" 
diagram, except of course for the extra thermal factor n(Ei). Thus the sum of all these diagrams, i.e., the diagrams 
that have only the i-th line snipped closed, adds up to n(Ei)Da(p, E). On the other hand, if the snipped line forms 
an external thermal leg pair (i.e., we have a "open" snipping, as in Figs.^d and^Jf), we again have an extra thermal 
factor n(Ei), but now the rest of the expression differs from that for the "un-snipped" graph in the sign of the energy 
Ei. This is so because, for an open snipping, the energy Ei moves from E cut to E inc , as can be gathered from rule 

^- 

Let x symbolize a variable, and let S x be the operator that acts on functions of x, changing the sign of the argument 
x, according to: 

S x f(x) : f(-x). 

In terms of the reflection operator S x , we can write the sum of all the time-ordered diagrams with only the i-th line 
snipped open as 

n(E l )S l D ( P , E), 

where we have written Si := Se { to avoid cluttering the notation. So the full contributions of the diagrams in which 
only the i-th line is snipped can be written as 



n(E i )(l+S i )D (p,E). 
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The analysis above can clearly be generalized to add up the contribution of the graphs with more than one snipped 
line. Taking into account that only connected graphs are allowed by the OFPT-rules (so that one is allowed to snip 
at most L internal lines, where L is the number of independent loops), we arrive at the following result: 

Theorem 1. The OFPT-rules admit the following mathematical representation: 

D R (p,E,T) = d(E,T)D (p,E), (55) 
where 0(E,T), the thermal operator, is given by 



6(E,T) :=l+Y / n(E l )(l + S l )+ ]T n(E n )n(E l2 )(l + S n )(l + S l2 ) 

i=l (iiM) 
L 

(ii, ...,»!,) 1=1 



(56) 






iP 






iP 





(a) 



(b) 





(c) 



(d) 





iP 










iP 






(e) 





iP 






/Ex 


iP 







(f) 



FIG. 4: An example of the diagrams which appear in the OPFT-rules. 
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Here the indices i\, 12, ■ ■ ■ run from 1 to I (the number of internal propagators) and the symbol . . . , ik) stands for 
an unordered k-tuple with no repeated indices. The primes on the summation symbols imply that we are to exclude 
from the sums those tuples (ii,. .. ,ik) such that if we snip all the corresponding lines .. ,ik then the graph becomes 
disconnected. 
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